Entanglement and decoherence in near-critical qubit chains 
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We study the problem of environmentally-induced decoherence in a near-critical one-dimensional 
system of N S> 1 coupled qubits. Using the Jordan- Wigner fermion representation of the qubit 
operators we identify the decoherence rates relevant for the two-qubit reduced density matrix. We 
find that a desirable onset of massive shared entanglement in the near-critical regime comes at the 
expense of decoherence which also tends to increase as the system is tuned towards criticality. Our 
results reveal rather contradictory general requirements that future designs of a qubit chain-based 
quantum information processor will need to satisfy. 



Besides establishing a number of other, previously un- 
explored, connections, the rapidly developing field of 
quantum computing brings about new links between sta- 
tistical mechanics of quantum spin systems and quantum 
information theory. 

This relationship was originally prompted by the idea 
of using electron or nuclear spins as natural candidates 
for individual physical qubits and constructing a practi- 
cal quantum register for implementing various quantum 
protocols. Moreover, it is believed to have an even greater 
potential, as far as practical control of quantum entan- 
glement and coherence is concerned. 

To this end, a number of authors have put forward 
the idea of constructing novel types of logical qubits 
which can enjoy an exceptionally high degree of coher- 
ence, thanks to the robust intrinsic correlations present 
in quantum states of many interacting spins. Instead 
of considering the qubit interactions as a nuisance to be 
rid of, these proposals seek to utilize them, focusing on 
topological and, therefore, intrinsically coherent Majo- 
rana fermion (edge) or anyon-like (bulk) excitations of 
some strongly correlated spin states proposed as a plat- 
form for implementing fault-tolerant computations with 
automatically built-in error correction codes. 

However, while offering a potentially ideal protection 
against environmentally induced decoherence, these pro- 
posals require a macroscopically large number of inter- 
acting physical qubits in order to construct only a few 
(topo)logical ones, and they may also face general diffi- 
culties with efficient information encoding and readout. 

Alternative, more memory-efficient, approaches to con- 
trolling decoherence, such as creation of decoherence-free 
subspaces (DFS) [1] by means of dynamical decoupling 
(" bang-bang" pulse control) [2] , often invoke certain sym- 
metries of the system-to-environment couplings and/or 
rely on the possibility of switching the qubits' interac- 
tions on and off at will, although, in practice, either con- 
dition may not necessarily be easy to meet. 

Regardless of the actual physical makeup of the qubits, 
most of the the previous analyses have been based on the 
generic spin-1/2 Hamiltonian 
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h= ]T £fl?(t)af+£J5(i)s?sy) (i) 

a=x,y,z i—1 k,j=1 

and an abstract quantum computing protocol has been 
thought of as a sequence of short pulses implement- 
ing one-qubit Bf{t) = B?8(t)9(T - t), and two-qubit, 
J?.(i) = J?.0(t)0(r - t), gate operations. 

It has been recognized, however, that any realistic 
qubits would interact with each other not only during the 
externally controlled gate operations but also during the 
idling periods due to their unwanted short-ranged (such 
as exchange, J?- cx exp(— const\i — j|)) or long-ranged 
(such as dipolar, Jf, cx — j| 3 ) static couplings. 

In this Letter, we study the possibility of utilizing these 
residual interactions for creating and maintaining quan- 
tum entanglement between the individual qubits that are 
also exposed to a noisy environment ("bath"). 

As one analytically tractable, yet sufficiently non- 
trivial, example we consider a system of N 1 
qubits with the topology of a one-dimensional (ID) 
chain and time-independent nearest-neighbor couplings, 
J°- = — J a Sj_i±i, while describing the bath by a fluc- 
tuating component of the local magnetic field, Bf(t) = 
B a + hf(t). 

Specifically, we focus on the two-parameter set of the 
Hamiltonians with B z = B, J x ^ = J(l ± 7), B x > y = 
J z — which is exactly solvable via the Jordan- Wigner 
(JW) transformation from spin-1/2 operators to spinless 
ID fermions 
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iS! = X il[e i ^ (2) 



Expressed in terms of the JW fermions, Eq.(l) takes on 
the form 

1 N 

H = - ^2 ( B xlxi - JXiXi+i - JjXiXi+i ) + h.c. (3) 

i=l 

that can be diagonalized by virtue of the Fourier transfor- 
mation Xj — ^2k elkaj Xk (here k = ±2irn/aN is the JW 
fermion momentum inversely proportional to the lattice 
spacing a), followed by the subsequent Bogoliubov ro- 
tation x±k = cosdkip±,k i sin0fc-0^ k through the angle 
defined by the relation tan 20^ = J7 sin ka/(B — J cos ka) 
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which turns (3) into a sum H = ^ ± k e k ip± k 4'±,k over 
single-fermion eigenstates ip+ lk and their charge conju- 
gates ("holes") ipl_ k with the dispersion 

e k = yj J 2 j 2 sin 2 ka+(B-J cos ka) 2 (4) 

At < 7 < 1 the Hamiltonian (3) belongs to the uni- 
versality class of the Ising model in transverse field, fea- 
turing a pscudo-relativistic low-energy dispersion e k ~ 
\Jv 2 k 2 + A 2 . The spectral gap (JW fermion "mass") 
A = \B — J | vanishes at the critical point A = B/J= f, 
whereas the velocity v = Jja remains finite. By contrast, 
at 7 = corresponding to the XY"-model the fermion dis- 
persion is quadratic at small momenta. 

In the thermodynamic (N — ► oo) limit, the model (3) 
demonstrates an onset of massive entanglement upon ap- 
proaching its critical point [3] , as quantified by the con- 
currence Cij = max[0, n — r 2 — r 3 — r&] defined in terms 
of the (numbered in the descending order) square roots 
r a of the eigenvalues of the product between the two- 
qubit reduced density matrix (RDM) pij = Tr-^p and 
its time-reversal conjugate of ® v V jP*j&\ ® a v - [4]. 

Near the critical point the nearest- and next- 
nearest-neighbor concurrences increase and their deriva- 
tives exhibit logarithmic singularities dCi^±\/d\ oc 
d 2 Ci } i±2 1 d\ 2 oc In | A — A c | with A c w 1 [3], consis- 
tent with the effective description in terms of the free 
(pseudo)relativistic ID fcrmions with mass A oc |A — A c |. 

The tendency of the correlated spin- liquid ("resonant 
valence bond") states to become more and more entan- 
gled in the critical (gapless) regime has to be distin- 
guished from the robust singlet correlations character- 
izing gapful "valence bond solid" states where, instead of 
being shared by many qubits, the entanglement is con- 
fined to their designated pairs. 

In light of the fact that distrubuted entanglement is 
considered an important resource for quantum informa- 
tion processing, one may be inclined to conclude that 
in order to take a full advantage of the increasing qubit 
correlations the system has to be tuned into the criti- 
cal (gapless) regime, e.g., by adjusting the uniform field 
{B —* J). 

However, any discussion of the potential benefits of 
the desired massive entanglement between the qubits 
would be incomplete without assessing the possibility of 
controlling their (in contrast to the former, unwanted) 
entanglement with the bath which constitutes a source 
of environmentally-induced decoherencc. In this regard, 
one of the key questions is whether or not this can be 
done by virtue of the same inter-qubit couplings that 
have facilitated the onset of the critical entanglement in 
the first place. 

In what follows, we choose the qubit-bath coupling 
Hamiltonian in the form X = Y^i=\ &i M (*)■ Unlike 
in the case of non-interacting qubits, the latter causes 



not only pure dephasing (phase errors) but also relax- 
ation (bit-flip errors), since [H, X] ^ 0. A spatial inho- 
mogeneity of a generic D-dimensional dissipative envi- 
ronment comprised of bosonic modes (photons, phonons, 
spin waves, etc.) with a dispersion U q oc requires one 
to use a dynamic spectral function 

N f oo 

D ab (uj,q)= V Im dte^*'^ < h^{t)h b j{0) > 

oc 5 ab Lu a+2 - D S(cu 2 - n 2 q ) (5) 

where the vector fy of length a\i — j\ is parallel to the 
chain, instead of a total spectral density J D zz (lj, q)d D q. 

Motivated by the results of Refs. [3] pertinent to the 
pairwisc entanglement and considering the special impor- 
tance of two-qubit gates in all the previously proposed 
quantum protocols, we will be primarily concerned with 
the two-qubit RDM which can be expanded over the one- 
and two-spin correlation functions 

Pij (t) = li ® lj + Y,(< Sf(t) >1 1 <E) a] + 

a 

< S?(t) > at ® 1,-) + < St(t)S b (0) > a? ® a) (6) 

a. b 

In particular, we will be interested in the off-diagonal 
RDM element 

<lVfa(t)\U>=<S+(t)S7(0) >= 

< Xi(t) n e -"^ (t)x ™ (t) J! e^(°>*»< o >xt(0) > (7) 
fe=l 1=1 

which, alongside the rest of the RDM, becomes 
translationally-invariant, <JT \f>ij(t)\ tl>= PtT;Ti(* — 
j, t), at qubit separations 1 <C \i — j\ «JV. 

Despite the fact that the dynamical correlator (7) ap- 
pears to be rather difficult to compute even in the noise- 
less case [5] , it gets affected by the noise in much the same 
way as the JW fermion propagator < Xi(*)xj(0) >; f° r 
the effect of the stochastic field h z (x, t) on the "charged" 
JW fermion operators \i an d x] dominates over that on 
their "neutral" products XmXm- 

In turn, the amplitude (7) can be expressed in terms of 
the (retarded) matrix-valued Green function G R (k 1 t) — 
9(t) < $(k,t)^(k,0) > expandable over the Pauli ma- 
trices fj acting in the Nambu (particle-hole) space of 
spinors "J/ = (ip+, tp-)- in the presence of a random field 
h z (x,t) and in the continuum limit, this Green function 
obeys the equation 

[id t -T 3 e k +fj k h z (x,t)]G R (k,t\h(x,t))=6(t) (8) 

where we introduced the matrix f] k = fs7] +tk + T\7]-^ k 
with the coefficients r] +ik = cos20fc and rj^^ k = sin26 k . 
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For a given configuration of the random field, an ap- 
proximate formal solution to Eq.(8) can be cast in the 
form (see Ref. [6] for the technical details of this tech- 
nique) 



|- J dse is ^- 



(9) 



to be interpreted as some kind of a "sub-exponential" 
(especially, temperature-independent) decoherence. 

In the non-interacting (J = 0) limit it is instructive to 
compare Eq.(10) to the well-known exact solution for the 
entire A-qubit density matrix [8] 



p(°)(i)=/3(0)exp 



1 2 J (2tt 



qdto 1 — cos tot 



D+l 



(12) 



exp [i 



The Fourier transform G R (x,t\h) ofEq.(9) represents the 
time evolution operator which incorporates the effect of 
" bremsstrahlung" (multiple emission/absorption of soft 
bosonic modes to/from the bath) and the associated re- 
coil of the JW fermions. 

With the operator of time evolution at hand, one 
can now construct the JW fcrmion density matrix 
pjw(x - y,t\h) = J dzdwG±±(x - z, t\h)p t j W (z - 
u>)G^_|_(w — y, — t\h) where the thermodynamic density 
matrix pf w {x) = J e lkx dk[l - tanh(e fc /2T)]/(47r) de- 
scribes the system of free JW fermions in equilibrium 
with the bath at temperature T. 

Upon performing a Gaussian statistical average of 
Eq.(9) over h z (x,t), one arrives at the following expres- 
sion for the JW fermion density matrix of a noisy near- 
critical qubit chain of length L = Na 



Pjw 



(x,t)= J L dyp%(y) [ ^e ik ^ cos 2 e k i 



2vr 



1 ^ f d D qduj 1 - cos(w - e k ± e k+q )t 



to 



T& tk D"(u,, g)(coth ^ - tanh ^^)(1 - co S q(x - y))} 



2T 



2T 



In the absence of the qubit-bath coupling, the kernel of 
the y-integration in Eq.(10) reads 



U(z,t) 



S(z) 



4tt S-j 



z±2vt 
z =F 2vt 



K^Ay/z* 



4vH 2 ) 
(11) 



where z = x — y+iO and K\(w) is the Macdonald function 
of the I s * kind. At A = and max[z, vt] < L it further 
reduces to the expression U(z,t) — [26 (z) + J2±3( z ^ 
2vt)]/4 which describes "ballistic" spreading of the ini- 
tial entanglement at the speed v away from the i th and 
j th qubits in both directions. 

We note, in passing, that, owing to its unitary na- 
ture, a spatial/temporal decay of the RDM of any ex- 
actly solvable noiseless qubit chain, such as Eq.(3) or the 
ferromagnetic Heisenberg chain studied in Ref. [7] , is not 
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where of and a i act upon the input (t = 0) density ma- 
trix p(0) on the left and from the right, respectively. 

By tracing out in Eq.(12) all but the i th and j th qubits, 
one readily obtains the formula for P[]-;-\i{i — j,t) which 
almost (see below) exactly coincides with Eq.(10) for 
v = 9 k = 0, thereby providing further support for link- 
ing the decoherence properties of the two-qubit RDM to 
those of the JW fermion density matrix pjw(x, t). 

In contrast to Eq.(10), however, in Eq.(12) the Fermi 
distribution function for the JW fermions is missing, con- 
sistent with the fact that in the non-interacting limit the 
qubit system lacks a well-defined temperature. As a re- 
sult, the frequency integral in Eq.(12) receives contribu- 
tions from the bosonic modes with energies to > T that 
give rise to the notorious T — decoherence, the latter 
being a ubiquitous feature of any Caldcira-Leggett-type 
(single-qubit) model which is oblivious to the Pauli ex- 
clusion principle. 

As follows from Eq.(12), under the conditions of 
"collective decoherence" (i.e., when the bath coherence 
length is much greater than L, hence e qXij s=s 1) the ma- 
trix element P^.^{i — j, t) demonstrates no decay at all, 
in agreement with the fact that the Hilbcrt space of the 
non-interacting qubits possesses a DFS spanned by the 
states that nullify the operator ~~ <?i) 

In fact, the validity of the above conclusion hinges on 
the possibility of characterising the bath solely in terms 
of the momentum-integrated spectral density. In con- 
trast to the case of non-interacting qubits, however, any 
non-flat JW fermion dispersion (e k ^= const) results in 
a non-trivial cf-dependcncc of the integrand in Eq.(10), 
thus making it difficult to positively establish the exis- 
tence of any DFS even under the most friendly collective 
decoherence conditions. 

Then one can readily see that for a generic encod- 
ing it is the slowest-decaying off-diagonal matrix ele- 
ment Pii;ii{i — j,t) that gets to control the long-time 
behavior of the concurrence Cij(t) and other entangle- 
ment quantifiers such as purity Pij(t) = Tr[p 2 j(t)] or 
fidelity Fij{t) = Tr[p t3 {t)e- lHt p l3 {0y Ht ]. 

Assuming that Eq.(10) retains its basic structure 

pjw{x,t) = [ dyU(x-y,t)p JW (y,0)e- r ^-y^ (13) 
Jo 
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for an arbitrary (not necessarily thermodynamic) input 
density matrix pjw(x,0), we now consider a practically 
important example of the isotropic D-dimensional bath 
composed of acoustic {j3 = 1) bosonic modes propagat- 
ing at a speed c >• v (albeit resulting in no significant 
loss of generality the latter condition does simplify the 
following analysis). 

For an input state that consists of closely spaced en- 
tangled qubit pairs (pjw(x,0) — > for x > a), Eq.(13) 
allows one to read off the proper dccoherence rate directly 
from the exponential attenuation factor (" influence func- 
tional") 

. . f°° , „ 1 - cos(w - A ± A)t 

r ( r )t) «y o ( ,_ A± A) 2 

, sin(wr/c) w , ui , u> — A N , „ 

Vl, k , (1 - ^^Mh - - tanh _) (14) 

where k* ~ 1/r. It is also worth noting that for a broadly 
distributed entanglement (/?jvi/(a;, 0) « const) the deco- 
herence rate can be more difficult to deduce. 

Representing the total decoherence rate (14) as a sum 
r = cos 2 26k,T+ + sin 2 26^ T_, we separate between the 
processes of scattering with almost no energy exchange 
(iv ~ 0), thus regarding T + as a loose analog of the "pure 
dcphasing" rate in the conventional singlc-qubit problem, 
and finite (lo w 2A) energy transfers characterized by the 
"relaxation" rate T_. 

Consistent with this interpretation, T_ also appears 
to control the (vanishing with decreasing J and/or in- 
creasing A, see below) decay rate of the diagonal ma- 
trix elements of Pij(t), as manifested by the correlator 
< S*(t)Sj(0) > related to the statistically averaged JW 

fermion amplitude < xl(*)Xi(*)xj(0)Xj(0) > 

Performing the integration in Eq.(14), we obtain 
r + (r, t) cx T(min[r/c,t]) 2 ~ a at long times/distances 
(min[r/c,t] > 1/T). In the opposite limit (max[r/c,t] < 
1/T) T + (r,t) behaves as cx T 3+a r 2 t 2 , while at interme- 
diate scales (min[r/c,t] < 1/T < max[r/c,t]) T + (r,t) cx 
T 1+a min[r/c,t] 2 . 

In the critical (A = 0) regime r_(r, t), too, is given 
by the above asymptotics, although a finite gap lim- 
its their applicability to max[r/c,t] < 1/A, while at 
min[r/c,t] > 1/A Eq.(14) yields T-(r,t) cx TtA a ~ 1 . 

Thus, at a = 1 corresponding to the Ohmic regime 
in the conventional spin-bath model the total decoher- 
ence rate F turns out to be essentially independent of 
A, as long as A < T, thus implying that at sufficiently 
high temperatures the quality of pairwise entanglement 
can be largely insensitive to the proximity to the critical 
point (the speed v of the ballistic entanglement spread- 
ing does decrease with increasing A, though). However, 
in the case of a sub-Ohmic (0 < a < 1) environment the 
decoherence rate (14) as a function of the gap develops a 
minimum at 



A- (l-af/^r 1 (15) 

Therefore, for a given time interval t between consecu- 
tive gate operations one can reduce the environmentally 
induced decoherence by tuning the system farther away 
from criticality and opening up the spectral gap (15). 

It also follows from Eq.(14) that, rather naturally, the 
decoherence rate decreases at smaller spatial separations 
between the qubits (large fc*), thus further emphasizing 
the special importance of the (next-) nearest-neighbor 
entanglement [3]. 

In summary, we studied the quality of time evolution 
between gate operations ("quantum memory") in a ID 
array of TV ^> 1 coupled qubits which are also subject to 
a spatially non-iniform (sub-)Ohmic dissipative environ- 
ment. We found that the shared quantum entanglement 
that reaches its maximal attainable value at the criti- 
cal points of the qubit-chain Hamiltonian [3] is generally 
accompanied by a concomitant increase in the entangle- 
ment between the qubits and the environment. There- 
fore, the requirement of preserving the qubits' entan- 
glement over a certain idling time between consecutive 
gates can be better fulfilled away from criticality and/or 
at smaller separations between the pair of qubits consti- 
tuting a gate. 

Viewed in the context of their potential applications 
to quantum information processing, our results suggest 
that, in contrast to the expectation expressed in the first 
of Refs. [3], the mere onset of the long-range correla- 
tions in the near-critical regime does not guarantee that 
an interacting multi-qubit system would autmatically be- 
come more robust against noise. In fact, the inter-qubit 
couplings appear to control both the onset of a critical 
behavior and the dependence of the decoherence rates on 
the proximity to a critical point. Therefore, when choos- 
ing the parameters of the qubit Hamiltonian (1), one will 
have to optimize between the rather contradictory crite- 
ria in order for a qubit-chain prototype of a quantum reg- 
ister to achieve its target performance in terms of both, 
the high degree of entanglement and reliable coherence 
control. 
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